Transformation seismology: composite soil lenses for steering surface elastic Rayleigh waves.

Metamaterials are artificially structured media that exibit properties beyond those usually encountered in nature. Typically they are developed for electromagnetic waves at millimetric down to nanometric scales, or for acoustics, at centimeter scales. By applying ideas from transformation optics we can steer Rayleigh-surface waves that are solutions of the vector Navier equations of elastodynamics. As a paradigm of the conformal geophysics that we are creating, we design a square arrangement of Luneburg lenses to reroute Rayleigh waves around a building with the dual aim of protection and minimizing the effect on the wavefront (cloaking). To show that this is practically realisable we deliberately choose to use material parameters readily available and this metalens consists of a composite soil structured with buried pillars made of softer material. The regular lattice of inclusions is homogenized to give an effective material with a radially varying velocity profile and hence varying the refractive index of the lens. We develop the theory and then use full 3D numerical simulations to conclusively demonstrate, at frequencies of seismological relevance 3–10 Hz, and for low-speed sedimentary soil (v(s): 300–500 m/s), that the vibration of a structure is reduced by up to 6 dB at its resonance frequency

Resonant metalenses for flexural waves in plates.

The dispersion curves of a cluster of closely spaced rods supported by a thin plate are characterised by subwavelength bandgaps and slow group velocities induced by local resonance effects. A recent analytical study [Williams, Roux, Rupin, and Kuperman (2015). Phys. Rev. B 91, 104307], has shown how the slow velocity branch depends, amongst other parameters, on the height of the rods that make up the cluster. Such metamaterial, offering easy-to-tune spatial velocity gradients, is a perfect candidate for building gradient index lenses such as Luneburg, Maxwell, and 90° rotating. Here theoretical results are combined with numerical simulations to design and test metalenses for flexural waves. The lenses are obtained by tuning the height of the cluster of rods such that they provide the required refractive index profile. Snapshots and videos from three-dimensional numerical simulations in a narrow band centered at ∼4 kHz are used to analyse the performances of three types of gradient index metalens (Luneburg, Maxwell, and 90° rotating)

Fluid-loaded metasurfaces

We consider wave propagation along fluid-loaded structures which take the form of an elastic plate augmented by an array of resonators forming a metasurface, that is, a surface structured with sub-wavelength resonators. Such surfaces have had considerable recent success for the control of wave propagation in electromagnetism and acoustics, by combining the vision of sub-wavelength wave manipulation, with the design, fabrication and size advantages associated with surface excitation. We explore one aspect of recent interest in this field: graded metasurfaces, but within the context of fluid-loaded structures. Graded metasurfaces allow for selective spatial frequency separation and are often referred to as exhibiting rainbow trapping. Experiments, and theory, have been developed for acoustic, electromagnetic, and even elastic, rainbow devices but this has not been approached for fluid-loaded structures that support surface waves coupled with the acoustic field in a bulk fluid. This surface wave, coupled with the fluid, can be used to create an additional effect by designing a metasurface to mode convert from surface to bulk waves. We demonstrate that sub-wavelength control is possible and that one can create both rainbow trapping and mode conversion phenomena for a fluid-loaded elastic plate model.Comment: 13 pages, 10 figure

Building Merger Trees from Cosmological N-body Simulations

Although a fair amount of work has been devoted to growing Monte-Carlo merger trees which resemble those built from an N-body simulation, comparatively little effort has been invested in quantifying the caveats one necessarily encounters when one extracts trees directly from such a simulation. To somewhat revert the tide, this paper seeks to provide its reader with a comprehensive study of the problems one faces when following this route. The first step to building merger histories of dark matter haloes and their subhaloes is to identify these structures in each of the time outputs (snapshots) produced by the simulation. Even though we discuss a particular implementation of such an algorithm (called AdaptaHOP) in this paper, we believe that our results do not depend on the exact details of the implementation but extend to most if not all (sub)structure finders. We then highlight different ways to build merger histories from AdaptaHOP haloes and subhaloes, contrasting their various advantages and drawbacks. We find that the best approach to (sub)halo merging histories is through an analysis that goes back and forth between identification and tree building rather than one which conducts a straightforward sequential treatment of these two steps. This is rooted in the complexity of the merging trees which have to depict an inherently dynamical process from the partial temporal information contained in the collection of instantaneous snapshots available from the N-body simulation.Comment: 19 pages, 28 figure

Experimental and Numerical Evidence of the Clustering Effect of Structures on Their Response during an Earthquake: A Case Study of Three Identical Towers in the City of Grenoble, France

In this article, interpretation of an equivalent to a macroseismic intensity survey, performed in three identical stand-alone buildings located in Grenoble, France, after an M L 4.1 earthquake, reveals a clustering effect, resulting in different levels of perception of seismic loading by inhabitants. The clustering effect is confirmed using numerical simulation; the variation of the seismic response of the building in the middle of the cluster depends on the azimuth of the seismic source relative to the building cluster. The major effect is the splitting of its resonance frequency, accompanied by a decrease in vibration amplitude. We conclude that clustering has an impact on urban effects, calling into question the validity of seismic design, which considers buildings in urban areas as stand-alone constructions, and the interpretation of macroseismic inten- sity surveys conducted in dense urban areas

Extension of tumor fingers: A comparison between an individual-cell based model and a measure theoretic approach

The invasive capability is fundamental in determining the malignancy of a solid tumor. In particular, tumor invasion fronts are characterized by different morphologies, which result both from cell-based processes (such as cell elasticity, adhesive properties and motility) and from subcellular molecular dynamics (such as growth factor internalization, ECM protein digestion and MMP secretion). Of particular relevance is the development of tumors with unstable fingered morphologies: they are in fact more aggressive and hard to be treated than smoother ones as, even if their invasive depth is limited, they are diffcult to be surgically removed. The phenomenon of malignant fingering has been reproduced with several mathematical approaches. In this respect, we here present a qualitative comparison between the results obtained by an individual cell-based model (an extended version of the cellular Potts model) and by a measure-based theoretic method. In particular, we show that in both cases a fundamental role in nger extension is played by intercellular adhesive forces and taxis-like migration

Adhesion and volume constraints via nonlocal interactions determine cell organisation and migration profiles

The description of the cell spatial pattern and characteristic distances is fundamental in a wide range of physio-pathological biological phenomena, from morphogenesis to cancer growth. Discrete particle models are widely used in this field, since they are focused on the cell-level of abstraction and are able to preserve the identity of single individuals reproducing their behavior. In particular, a fundamental role in determining the usefulness and the realism of a particle mathematical approach is played by the choice of the intercellular pairwise interaction kernel and by the estimate of its parameters. The aim of the paper is to demonstrate how the concept of H-stability, deriving from statistical mechanics, can have important implications in this respect. For any given interaction kernel, it in fact allows to a priori predict the regions of the free parameter space that result in stable configurations of the system characterized by a finite and strictly positive minimal interparticle distance, which is fundamental when dealing with biological phenomena. The proposed analytical arguments are indeed able to restrict the range of possible variations of selected model coefficients, whose exact estimate however requires further investigations (e.g., fitting with empirical data), as illustrated in this paper by series of representative simulations dealing with cell colony reorganization, sorting phenomena and zebrafish embryonic development

Collective migration and patterning during early development of zebrafish posterior lateral line

The morphogenesis of zebrafish posterior lateral line (PLL) is a good predictive model largely used in biology to study cell coordinated reorganization and collective migration regulating pathologies and human embryonic processes. PLL development involves the formation of a placode formed by epithelial cells with mesenchymal characteristics which migrates within the animal myoseptum while cyclically assembling and depositing rosette-like clusters (progenitors of neuromast structures). The overall process mainly relies on the activity of specific diffusive chemicals, which trigger collective directional migration and patterning. Cell proliferation and cascade of phenotypic transitions play a fundamental role as well. The investigation on the mechanisms regulating such a complex morphogenesis has become a research topic, in the last decades, also for the mathematical community. In this respect, we present a multiscale hybrid model integrating a discrete approach for the cellular level and a continuous description for the molecular scale. The resulting numerical simulations are then able to reproduce both the evolution of wild-type (i.e. normal) embryos and the pathological behaviour resulting form experimental manipulations involving laser ablation. A qualitative analysis of the dependence of these model outcomes from cell-cell mutual interactions, cell chemical sensitivity and internalization rates is included. The aim is first to validate the model, as well as the estimated parameter values, and then to predict what happens in situations not tested yet experimentally. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'